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1 Introduction 

As an important application of chaotic dynamical systems, chaos-based secure commu- 
nication and cryptography attracted continuous interest over the last decade. It studies 
methods of controlling deterministic systems with chaotic behavior. Moreover, it is easy 
to notice the possibility of substantial variation of the characteristics of chaotic systems 
by relatively small variations of their parameters and external actions. A method of 
transmitting information using chaotic signal was proposed by A.S.Pikovsky and M.I. 



Rabinovich ( jl2l ) using the differential system 



xi = —lyiXi + hx2 + X2X3 
X2 = hxi — 1/2X2 — X1X3 
X3 = -1^3X3 + X1X2. 

In this paper we will consider the Rabinovich system: 

Xi = X2X3 

±2 = -X1X3 , (1.1) 
X3 = X1X2 

and we will analyze some global properties, the local study of stationary points, compatible 
Poisson structures and corresponding tri-Hamiltonian systems are also discussed. 

A Hamiltonian equation is called tri-Hamiltonian if it admits two Hamiltonian repre- 
sentations with compatible Poisson structures 

dx ~ ~ _ _ 

— = J\/H = JVH = JVH, 
dt 

where J, J and J are three Hamiltonian matrices (of the form [{xjjXj}] and {■, ■} is the 
Poisson structure) and they are also compatible. 
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The Hamiltonian formulation is important in mathematics, physics and also in other 
branches of natural science. From another point of view, we considered the analysis of the 
revised dynamical system, the analysis of the dynamical system with distributive delay 
variables and the analysis of the fractional dynamical system. 

Revised dynamical system, with distributive delay, associated to system fll.ip . allow 
the description of new crypting methods. 



2 The analysis of classical Rabinovich differential 
equations 



2.1 Geometrical properties of the system (11.11) 

In this subsection we will present some dynamical and geometrical properties, from 
geometrical mechanical point of view, ((0),([3))- 

Proposition 2.1 The dynamics ( li. ij) have the following Hamilton-Poisson realizations: 
(i) {M:\P\hi), i = 1,2,2 where 



/ii(xi,a;2,X3) = -{x\ + xl)] 
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/i2(xi,a;2,X3) = xl + xl; 

h^{Xi,X2,X^\ 



(it) (M3,Pif3^^a), where 



PuV = <yPi+PP2 + lP3 





—ax3 



axs 




{a-l + ^)x2 (f-i)xi 



(a - I + |)x2 

-(f-i)^i 




and ha{xi,X2, x^,) = ^{xf + x^), for each a, /?, 7 G M, a 7^ 0. □ 

From direct computations, using the algebraic technique of Bermejo and Fairen 
we get the following results. 

Proposition 2.2 There exists only one functionally independent Casimir of our Poisson 
configurations: 



(]R^Pl) given by: ci(a;i, X2, X3) 



2 (^2 + ^3)5 
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(M^, P2) given by: C2{xi,X2, X3) = xl + xj; 
(]R^P3) given by: C3(xi, X2, X3) = xj + xj; 

a^ven by: c„^-y(xi, Xa, Xs) = -^(f - + ^1 



□ 



2.2 Stability problem 



From the analysis of the stationary points, using fjl3l ). we get the following statements 



Proposition 2.3 The stationary points e^{m,0,0), 62^(0, m,0) and e^(0,0,m), m G M 
have the following behavior: 

(i) e™(m, 0,0), m G M are spectrally stable; 

(a) e^(0,m, 0), m G M are unstable; 

(Hi) e™(0,0,m), m G M are spectrally stable.^ 

Proposition 2.4 The stationary points e™(m, 0,0) and e^(0,0,m) are nonlinear 
stable.n 

Now we will point out periodic orbits and heteroclinic orbits associated to Rabinovich 
system. When studying the existence of periodic solution, we will use (jsj). 

Proposition 2.5 The dynamics reduced to the coadjoint orbit (xi)^ + (^2)^ = rri^ has 
near the stationary points e™(m,0,0) m G M* at least one periodic solution whose period 
is close to T^.D 

\m\ 

Proposition 2.6 There exists four heteroclinic orbits between the stationary points 
e^(0, m, 0) and e^™'(0, — m, 0), m G M, m 7^ given by: 

Xi{t) = ±msech{mt) 

X2{t) = ±mtgh{mt) (2.1) 
X3{t) = ±m sech{mt) . 

These orbits belong to the planes X3 = ixi.D 

3 The metriplectic structure associated to 
Rabinovich system 

A Leibniz structure on a smooth manifold M is defined by a tensor field P of type (2,0), 
. The tensor field P and a smooth function h on M, called a Hamiltonian function, 
define a vector field Xh which generates a differential system, called Leibniz system. If 



3 



P is skew-symmetric then we have an almost simplectic structure and if P is symmetric 
then we have an almost metric structure. 

Let P be a skew-symmetric tensor field in M.^ of type (2,0), g a 2-symmetric tensor 
field and h G C°°(]R^). If P is a Poisson tensor field and g is a nondegenerate tensor field, 
then (]R'^,P, (/) is called a metriplectic manifold of the first kind ((jy), (jiol ). (jlll)). The 
differential system is given by 



Xi 



p 



j=l J j=l 

The Qij is compatible with h, and is given by: 

on. on 

9ij 



1,2,3. 



9v 



E 



dxk dxk ' 



dh dh 



z,j = 1,2,3. 



(3.1) 



(3.2) 



If P is a (almost) Poisson differential system on with Hamiltonian function hi and 



a Casimir function /i2, there exists a tensor field g such that (M^,P, (/) is a metriplectic 
manifold of second kind. The differential system associated with it is given by: 



Xi 



dxi 



1,2,3, 



(3.3) 



where 



Let 



S pc 



a 



- E 



dhi dh2 
dxk dxk ' 



9ij 



dhi dh2 
dxi dxj 



z,j = 1,2,3. 



(3.4) 



1,2,3 realizations of Rabinovich system of differential equations. 



with Hamiltonian functions ha, a = 1, 2, 3 and Casimir functions Ca, a = 1,2, 3, where: 



X3 
-X3 
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2\ — 2 _|_ 2\ 
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\X2 ■ 





-\X2 
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/ 2 I 2\ / 2 I 2\ 

-[X1 + X2), C2 = -[X2 + x^); 



h^ 



{xf + xD, C3 



(^2 + ^3) 



Using (I3IID, (1321), (B> and (I33D we get the following resuhs. 



(3.6) 



(3.7) 
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Proposition 3.1 (a) The metriplectic realization of the first kind of 5\) is given by 
(M^, Pi,gi) where: 

9n = -(a;2)^ 9I2 = -{xif, ^33 = 0; 
gu = X1X2, g\i = X1X2, gls = gli = 0, ^23 = 9I2 = 0- 
(b) The associated differential system is given by: 



(3i 



Xi = X2X3 + XiX2{xi - X2) 
X2 = -X1X3 + (Xi)^ 
X3 = X1X2. 

(c) The differential system / IJ.<^) has the following stationary points: 

e^KCO), e™(0,m,0), e™(0,0,m). 

(d) The matrix of the linear part of the system ( I J. i^j) in e™(m, 0,0), e^(0,m, 0), resp 
in e^(0,0,m) is given by: 
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— m 
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(e) The characteristic equation of Ai for in e^(m,0,0) is: 

Xi-X^ + m^{m + l)) = 

and so, we have two cases: 

(i) if m > —1, then e^{m,0,0) are unstable; 

(ii) if m < —1, then we have a limit cycle. 

(f) The characteristic equation of A2 for ^3. <S)) in e^(0,m,0) is: 

-A(A^ + m^A - m^) = 

and so, it can be easily seen that 62(0, m,0) are unstable. 

(g) The characteristic equation of A3 for ^3. <SI) in e|(0,0,m) is: 

A(A2 + m2) = 0. 

(h) In a neighborhood o/e^(0,0,m), m > there exists a limit cycled 

Proposition 3.2 (a) The metriplectic realization of the second kind of 1^3. 5\} is given by 
(M^, Pi,gi) where: 

gli = -(a;2)^ ^22 = 0, ^33 = 0; 
gl2 = X1X2, gli = 0, gl^ = X1X3, gl-^ = 0, ^23 = ^2X3, ^32 = 0- 
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(b ) The associated differential system is given by: 



Xl = X2X3 + Xi{{x2f + (x-sf) 

X2 = -xixs + X2X3 (3.9) 

X's = X1X2. 

(a) The differential system y\) has the following stationary points: 

er(m,0,0), e™(0,m,0), e™(0,0,m). 

(d) The matrix of the linear part of the system Ii3. 8\) in e™(m, 0,0), e^(0,m, 0), resp 
e^(0,0,m) is given by: 





'0 





" 
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m 





m 
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— m 


, A2 = 
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, resp A3 = 


—m 


m 










m 
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(e) The characteristic equation of Ai for \3.9^) in e^(m,0,0) is: 

A(A=' + m2) = 0. 

(f) The characteristic equation of A2 for \3.9^) in e^(0,m, 0) is: 

A(A2 - Am^ - m^) = 

and so, it can be easily seen that e^(0,m,0) are unstable. 

(g) The characteristic equation of A^ for ^3. 9^) in e™(0,0,m) is: 



\{\^ - \{m + m^) + m^) = 0. 



□ 



Remark 3.3 In an analogous way we can discuss the metriplectic realization of first kind 
of ISJj and (M. 



4 The differential systems with distributed delay 

Let us consider the product M'^ x M'^ = {(x, x) | x G M^, x G M'^} and the canonical 
projections TTj : x — , i = 1, 2. A vector field X G X(R^ x R^), satisfying the 
condition X{nlf) = 0, for any / G C°°(M^), is given by: 

" d 

X{S:,x) = J2Mx,x) — . (4.1) 

i=l * 

The differential system associated to X is given by: 

Xi{t) = X{i,x), i = l,2,3. (4.2) 
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A differential system with distributed delay, see (|2|) is a differential system associated 
to a vector field X e x M?) for which X{Tilf) = 0, V/ G C°°(M3), and it is given 

by (14.11) where x{t) is: 

x{t) = / k{s)x{t - s)ds (4.3) 
Jo 

k(s) is a distribution density. In the following we will consider the following densities: 
1. uniform: 



0, < s < a 

-, a < s < a + T 

0, s > a + T. 



(4.4) 



where a > 0, r > are fixed numbers. 
2. exponential: 



ka{s) = ae "^ 



3. Erlang: 

4. Dirac: 



kn,(s) = a se 



2 „„-as 



a > 0; 



a > 0; 



ka{s) =S{s-t), t>0; 



(4.5) 
(4.6) 
(4.7) 



The differential equations with distributed delay for Rabinovich system are generated 
by an antisymmetric tensor field P on x that satisfies the following relations: 

^(vTiVi, </2) = 0, P(7r;A, 7r*/2) = 

for all /i,/2 e C°°(M3). 

The differential equation with distributed delay is given by: 



x{t) = P{x{t),x{t))V,h{x{t),x{t)), 



(4. 



where x{t) = k{s)x(t — s)ds, and h E C 
Let 
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We define 



P{x,x) = '^^eiPi, with > 0, ^^£, = 1. 



2=0 



i=0 



(4.9) 



Let 



ho{x,x) = ^((xi)^ + (xa)^), = XiXi + ^(a:2)^ 

h2{x,x) = ^{xiY + X2X2,h3{x,x) = XiXi + X2X2. 



We define 

3 3 
h{x,x) = ^Sihi, with 5i > 0, ^5i = l. (4.10) 

1=0 i=0 

The Rabinovich differential equation with distributed delay is given by (14.81) with P 
and h given above by (14. 9 p and (I4.10p with initial value x{s) = (f){s), s G (— oo,0] where 
(j) : (-00,0] — > M^, G C°°(M=^). 

In what follows we consider the functions / G C°^(M'^ x M.^) given by: 

lo{x,x) = ^{{X2Y + {X3f),li{x,x) = X2X2 + ^{X3Y, 
hix.x) = l-{x2y + X3X3,l3{x,x) = X2X2 + X3X3. 



We define 



/(z, z) = ^ Sj/jj, with £4 > 0, ^£j = l. (4.11) 

'i=0 1=0 

Proposition 4.1 • The junction l{x,x) given by Ili4-ll\ ) satisfies the following rela- 
tion: 

VJ{x, x)P{x, x)V^f{x, x)=0, / G C°°{R^ X R3); (4.12) 

• The revised differential equations with distributed delay satisfies the following rela- 
tion: 

x{t) = P{x, x)'Vxh{x, x) + g{x, x)Vxl{x, x) (4-13) 
where x{t) = k{s)x{t — s)ds and g{x, x) is a 2-tensor field given by: 



g{x,x) = {gij{x,x)), 

/_ X dh(x,x) dh(x,x) . , . 

^-^"'"^ = d^^'^' 

^ (dh{x,x)\^ 
g.,{x,x) = - ^ [-^^) ■ 



k=i,k^i 



□ 
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The revised Rabinovich system with distributed delay has the following form: 

Xi = {aiX-i + 04X3) (/9iX2 + (3^X2) + {P2X1 + P3Xi){f3iX2 + P4X2)a2X3 
±2 = -{aiXs + a4X3){P2Xl + (^3X1) - {P2X1 + P3Xi)a3X3 

X3 = 02X2 + a3X2){l32Xi + iSsXi) - {P2X1 + /33Xi)q!4X3 - (/3iX2 + /94£2)^a3a;3 

(4.14) 

where we considered the following notations: 

ai = Eq + El, a2 = Eo + £2, as = Ci + £2, = Ei + E3, = E3 + £2 

Pi = So + 61, ^2 = ^0 + ^2, /33 = Sl + S3, /34 = 62 + 63. 

Remark 4.2 The analysis of stationary points of the system Ili4-14\ ) is quite difficult, that 
is why we will present the main results for fractional Rabinovich differential system. 



5 Fractional Rabinovich differential systems 

Generally speaking, there are three mostly used definitions for fractional derivatives, 
i.e. Griinwald-Latnikov fractional derivatives, Riemann-Liouville fractional derivatives 
and Caputo's fractional derivatives, (([l|),(j9|)). Here we discuss Caputo derivative: 

— j x{t), a>0 (5.1) 

where m — l<a<m, m>l, (^j =^o...o^, /^is the Z?*'^ order Riemann-Lioville 
integral operator, which is expressed in the following manner: 

I^xit) = — — / (t - sf-^x{s)ds, P>0. (5.2) 



r(/5) Jo 

In this paper we consider that a G (0, 1). 

A fractional system of differential equations with distributed delay in M.^ is given by: 

D^x{t) = X{x{t),S:{t)), ae(0,l) (5.3) 

where x{t) = {xi(t),X2{t),X3(t)) G M^. 

The matrix associated to the linear part of the system (15.31) in the stationary point xq 
is given by the linear fractional differential system: 

D^u{t) = Au{t) + Bu{t), (5.4) 



where A = ) and B = (^|f ) 



x=xo 



The characteristic equation of (15.41) is: 

A{X) = det{X'' -A-k\\)B) (5.5) 
where fc^(A) = k{s)e~~^^ds and k is given by (14. 4p - ( 13171) . 
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Proposition 5.1 (0j 

1. If all the roots of the characteristic equation A (A) = have negative real parts, then 
the stationary point xq of Ii5.5\) is asymptotically stable. 

2. If k(s) is Dirac distribution, the characteristic equation is given by: 

A(A) = det(A"-A-e~^"5) = 0. (5.6) 

If T = 0, a ^ (0, 1) and all the roots of the equation det{XI — A — B) = satisfies 
I arg{X) |> then the stationary point xq is asymptotically stable. 

3. If a & (0.5, 1) and the equation det{XI — A — e~^'^B) = has no purely imaginary 
roots for any r > 0, then the stationary point is asymptotically stable. 

□ 

Let us consider a fractional 2-tensor field e X"(M=') x X'^{R^) and rf"/, d'^g G T){R^). 
The bilinear map 

[-, : C°°(M3) X C°°(M3) — , C°°(M3) defined by: 

[f,gr = B'^id''f,d^g), f,g eC^iR') 

is called the fractional Leibniz bracket. 

If is skew-symmetric, we say that (M^, [■, ■]") is a fractional almost Poisson manifold. 
For h G C°°(M^) the fractional almost Poisson dynamical system is given by: 

3 

Dfx,{t) = [x,{t),h{t)r, [x„ hr = J2 ^"^"r (5-7) 

Let be a skew-symmetric fractional 2-tensor field and a symmetric fractional 
2-tensor field g"' on M^. We define the bracket 
[-, : C°°(M3) X C°°(M3) , C°°(M3) by: 

[/, gr = P"(rf"/, d'^h) + g^idy, d-h), f,he C°°(M3). 

The 4-tuple {M.^,P",g°', [■,■]") is called fractional almost metric manifold. The frac- 
tional dynamical system associated to /i G C°°(M^) and 

3 3 

D^x,{t) = [x,{t),h{t)r, [x„ hr = + E ^r^D:. (5.8) 

i,j=l i,j=l 



Proposition 5.2 1. The fractional dynamical system ( [5. 7\ ) is given by: 

D^xiit) = X2{t)x^{t) 

Dfx2{t) = -xi(t)x3(t) (5.9) 
D^xsit) = Xi{t)x2{t) 
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2. The fractional dynamical system Ii5.8\) is given by: 



' D^x,{t) = ((aiX3(t) + a^X3it)){f3iX2it) + (3^i2{t)) 

+ l(32Xi{t) + (33Xi{t)){^iX2{t) + (34X2{t))c^2X3{t) 

< Dfx2{t) = -{aixsit) + a4X3{t)){(32Xi{t) + p3Xi{t)) - {P2Xi{t) + /^a^iW) 032:3 (t) 

Dfxsit) = {a2X2{t) + a3X2{t)){p2Xi{t) + p3Xl{t)) 
, -{p2Xl{t) + (33Xi{t))a4X3{t) - {(3iX2{t) + PiX2{t))''a3X3{t). 

(5.10) 

3. The fractional dynamical systems 1^5.91) and ( (5.101) have the stationary points 
e^(m, 0, 0), e^(0, m, 0) and e^(0, 0, m), m G M. 

4- The characteristic equations for Ii5.y\) are given by: 

• e^(m, 0, 0) : A"(-A2" + m^{m + 1)) = 0; 

• e^(0, m, 0) : A°(A2" + m^A^ - m^) = 0; 

• e^(0, 0, m) : A"(A2" + m^) = 0. 

5. r/ie characteristic equations for ( (5. iOI) are (72^ en 

• e5"(m, 0, 0) : A°(A2" + aA" + 6e~^"^ + c) = 0, a, fe, c G M; 

• e^(0, m, 0) : A"(A2" + oiA" + 6ie~^"^ + Cie-^^"^) = 0, ai, 61, Ci G M; 

• e^(0, 0, m) : A°(A2" + osA" + bae"^"" + Cse-^^"" + 4) = 0, 02, &2, C2, ^2 G M. 

In the second section, "The analysis of classical Rabinovich differential equations", 
we worked with a single-step method, Runge-Kutta, that means that we used only the 
information regarding the previous point for computing the successive point. For our 
illustrations we develop the Adams-Bashforth-Moulton predictor-corrector method. 

We can integrate numerically the set of differential equations fl5.9p with the Adams- 
Bashforth-Moulton method. To do so, we consider the relations we used for our simulation 
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xip{j + l) = 2:1(0) + 

X2{j + l) = 2:2(0) + 

X2p{j + 1) = 2:2(0) - 

2:3(^ + 1) = X3(0) + 



r(a 
1 



r(a 
1 



part, in Moulton method: 

1 ^ 

XiU + 1) = 2:1(0) + f^C^(^ik,j + I)x2{k)x3{k) + a{j + 1, j + l)x2p{j + l)x3p(j + 1)) 

fc=0 

j 

{J2KkJ + l)x2{k)xs{k)) 

k=0 

j 

(^ -a{k,j + l)xi{k)xs{k) - a(j + 1, j + l)xip(j + l)2:3p(j + 1)) 

A;=0 
j 

{Y,Kk,J + l)xi{k)x3{k)) 

k=0 

i 

{^a{kj + l)xi{k)x2ik) + a{j + l,j + l)xip{j + l)a:2p(i + 1)) 

fc=0 

j 

(^6(A;,j + l)xi(A:)x2(fc)), 



r(a 
1 



r(a 
1 



X3p(j + 1) = X3(0) + 



1 



r(a 



fc=0 



where 
and 



b{t, j + = 

a(i + l,j + l) = h 

a{OJ + l) = /i" 

aU + l,j + l) -- 



i = 0,l,...,n and j = 0, l,...,m 

„(j-^ + l)"-(j-z)" 



, (j - ^ + 1)"+^ + (j - ^ - 1)°+^ - 2(j - z) 
a{a + 1) 

r+'-(j-«)(j + i)" 



a+l 



a{a + 1) 



a(a + 1) 



We consider a graphic representation of Moulton method, for our system (11.11) . for the 
following two cases: 

(1) xi(0) = 0.001, X2(0) = 0.001, X3(0) = 6, a = 0.8 (Figure 1); 

(2) xi(0) = 0.001, X2(0) = 0.001, X3(0) = 6, a = 1 (Figure 2) 
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Figure 1 alpha=0.8 




6 Conclusions 



Until now we have an approach and also some solutions for metriplectic manifolds 
of the first and second kind, and also differential systems with distributed delay(for 
the first case presented here), and Rabinovich fractional differential system, with Dirac 
distribution(x(t) = x{t — r)).What we want to continue is to apply all other distributions 
for our three cases. 

Acknowledgements : The authors were partially supported by the Grant GNOSIS 
95GR 2007/2008. 
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